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Abstract 

The c — > w~f transition does not occur at the tree level in the Standard Model and is 
strongly GIM-suppressed at one loop. The leading QCD logarithms are known to enhance 
the amplitude by more than an order of magnitude. We point out that the amplitude in- 
creases further by two orders of magnitude after including the complete 2-loop QCD contri- 
butions. Nevertheless, AS = radiative decays of charmed hadrons remain dominated by the 
c — ► ddwy and c — ► sswy quark subprocesses. 
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1. Inclusive heavy flavor decay can be systematically analyzed with help of an expansion 
in inverse powers of the heavy quark mass Q. At the leading order in such an expansion, 
the inclusive hadronic decay rate is given by the perturbatively calculable free quark decay 
rate. For charmed hadron decays, this procedure can certainly be questioned since the charm 
quark mass is not really much larger than the QCD scale A. Nevertheless, some properties 
of charmed hadrons have been analyzed succesfully with help of HQET [[J. Therefore, it is 
also of interest to calculate inclusive decay rates of these particles at the leading order in the 
heavy mass expansion. 

In the present letter, we consider weak radiative charm decay |J. For definiteness, 
we restrict ourselves to AS" = processes. At the leading order in electroweak and strong 
interactions, there are three contributions to the AS = charm quark radiative decay: 
c — > wy , c — > udd'j and c — > uss'y. At higher orders in the strong interactions, there 
are virtual corrections to these decays as well as new decay modes with more gluons and 
quark- ant iquark pairs in the final state. 

In the absence of QCD, the c — > uy amplitude is enormously suppressed by GIM can- 
cellations and by small Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The leading 
logarithmic contribution turns out to be considerably larger, but is still afflicted with a small 
CKM coefficient. This raises a hope that higher order terms can be the most important ones, 
if they are not suppressed by small CKM factors. 

In this paper, we show that the c — > wy process is completely dominated by a two- 
loop term which has not been considered so far in perturbative analyses of charm decays. 
Furthermore, we point out that decays with more gluons which are also of higher order in 
the strong interactions have comparable amplitudes. 

The c — > wy channel cannot be separated kinematically from the accompanying modes 
by making a cut on the photon energy spectrum, because all the decay products are light. 
This is contrary to the b —>■ s'j decay which can be separated from the b — > scuy and other 
similar channels by selecting photon energies above the charm production threshold. In the 
c — > wy case, we might still hope to be able to separate this channel by looking at the shape 
of photon energy spectrum. However, the possibility of doing so depends crucially on how 
large the c — > wy amplitude is. 
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2. We begin with the c — > wy process and briefly review the known contributions which come 
at the (formally) lower orders of perturbation theory. All momentum invariants in a two- 
body decay are expressible in terms of masses. Therefore, we can represent the amplitude 
by a tree diagram with a single local vertex. For the c — > wy decay in the Standard Model, 
the relevant local interaction reads 



J int 



AG p 

A — — - m c (ua^P R c)F^, 



v/2 " 16?r 2 

where P R = (l+7 5 )/2 and we have neglected the u-quark mass. The coefficient A nontrivially 
depends on light quark masses (see below). Thus, it should not be misinterpreted as a Wilson 
coefficient. 
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Figure 1: One- loop diagrams generating the c — > wy transition. 



As any flavor changing neutral current process, the c — > wy amplitude arises in the 
Standard Model only at the one-loop level. The relevant diagrams are shown in fig. |l] and 
give the following contribution to the coefficient A: 



loop 
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q=d,s,i 



M. 
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The CKM factors in the above equation have very different orders of magnitude 



\v; d v ud \ ~ \v; s v us \ ~ 0.22 and \v; b v ub \ ~ (1.3 ± 0.4) 



x 10" 



(3) 



Consequently, |AA 1!oop | ~ 2 x 10 -7 . The extraordinary smallness of this number is due to 
the tiny factors {m q /My/) 2 for the light quarks and to the small CKM angles in the 6-quark 
contribution. 

Since the important suppression factors are independent of gauge couplings, it is possi- 
ble that higher orders in perturbation theory give dominant contributions to the radiative 
amplitude considered because they may not suffer the same dramatic supressions and are 
reduced only by powers of the gauge couplings. 

A natural first attempt to include higher orders is to resum short distance QCD corrections 
in the leading-logarithmic approximation, by analogy to the b — > sj decay where they 



bring sizeable enhancements pi. In order to systematically include these contributions, one 
introduces two effective hamiltonians 



H eff (M w > /i > m b ) 
H eff (m b > fi> m c ) 



AG 

~i E v^cmoi + cmo*] 

q=d,s,b 
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V 7 ^ q t^, s 



E V*V uq [Cx{n)0\ + CMOl + E CiMOi 



(4) 
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where 
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{u a -f IM P L q f 3)(q a J ,1 PLCf3), 
m c (u a a^P R c a )F^ v . 



q 



d, s, b 



d, s, b 



(6) 
(7) 
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The remaining operators Oj are given explicitly in eqn. (18) of ref. ||. We do not present 
them here because they are inessential for our discussion. 

As indicated in eqns. (^) and (||), each of the hamiltonians is valid in a different range of 
the renormalization scale fi. Now, we neglect all terms which are suppressed by additional 
powers of l/M w such as the ones in eqn. (|2|). Due to CKM unitarity, the operators 3 ,...,0 8 
are not generated by QCD renormalization at scales fi > mi,. 

Finding the coefficients Ci(fi = Mw) and performing the renormalization group evolution 
from fi = Mw to fi = m c is by now standard ||. Analogously to the b — > 57 analysis of 
ref. 0, we introduce an "effective" anomalous dimension matrix. It reads 
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For the renormalization group evolution from fi = m b to fi = m c in the c — > W7 case, one 
needs to substitute / = 4, Qi = Qd = — 5, Q2 = Q« = § and Q = 2Q U + 2Qd = §. 

For the evolution from /x = to fi = m b in the 6 — > 37 case, one would need to 
substitute / = 5, Qi = <5«, Q2 = Qd and Q = 2Q U + 3Qd 



|. In this case, the matrix 
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(0)eff 



given in Appendix A of ref. || would be reproduced. 



The resulting leading logarithmic contribution to A in eqn. (JI]) reads 

AA LLA = -V; b V ub C 7 e// (m c ), (9) 

where 

c e 7 ff (m c ) = f; (^t) 1 WW + *)] • ( 10 ) 

The coefficients Ci^ijiib) are given by 



_6_ , .12 

" 23 



2 \ a s (m fe ) / 2 \ cn s (m&) / 
and the numbers Xi, yi and Zj read 

= ( -fli' i§T> §» °' °- 6524 ' -0-0532, -0.0034, -0.0084), 
= ( -f|§ff, f§, f, 0, 0.8461, 0.0444, 0.0068, -0.0059), (12) 
*< = ( M> M> i' "§> 0- 3469 ' -0-4201, -0.8451, 0.1317). 

Note that the r.h.s. of eqn. ( |T0"D vanishes in the formal limit a s (m c ) — > a s (mi,), as it should. 

Taking m b = 5 GeV, m c = 1.5 GeV and a s (M z ) = 0.12 (which implies a s (m b ) / a s {m c ) ~ 
0.67), we obtain 

\AA LLA \ = [O.OOldK) + 0.055C 2 (m 6 )] \V^V uh \ = 0.060 \V: b V uh \ ~ (8 ± 3) x lO" 6 . (13) 

This result is more than an order of magnitude larger than the (formally) leading order one in 
eqn. (0): Including logarithmic QCD contributions replaces the powerlike GIM suppression 
factors (m q /Mw) 2 by a logarithmic function of m b /m c . 

This surprising enhancement was noted some time ago and discussed in detail in ref. 
Pf. However, our coefficient Cy (m c ) in eqn. flIU|) is 5 times smaller than the one obtained 



there. This difference arises mainly because the operators 0\ and 0\ have not been taken 
into account in ref. |J. Their presence above \i = m b is essential for cancellation of additive 
logarithmic QCD contributions at these renormalization scales. 

3. Although the leading logarithmic result (p~3|) is much larger than the purely electroweak 
contribution (0), there remains a strong cancellation between s and d loops whose CKM 
factors are very similar in magnitude but have opposite signs. Consequently, a suppression 



by V* b V ub is still present in eqn. fll3 



This cancellation is circumvented when the functional dependence of the amplitude on the 
s- and <i-quark masses becomes substantial. It turns out that this happens in the two-loop 
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Figure 2: Diagrams contributing to the two-loop c — > wy matrix element of 0\. 



matrix element of the effective hamiltonian in eqn. (|5|). It is given by the diagrams in fig. ^] 
with 2 and 2 operator insertions.0 

After using unitarity of the CKM matrix, we obtain the following contribution to the 
coefficient A: 



A = r;y us ^4^ C 2 {m c ) {f[(m s /m c ) 2 ] - f[(m d /m c ) 2 ]} + 0(V^V ub ) 



4tt 



(14) 



where 



c 2 {m ) = - (^AMeI 

2 I a s (m b ) 



23 / a s (m b ) \ 25 



1 faJMv) 

2 1 a s (m b ) 



12 

" 23 



12 

" 25 



(15) 



\a s (m c ) J 2\a s {m b )) \a s (m c ) J 

The function / can be extracted from an analogous computation performed recently for the 
6-quark decay M. It reads 

— |576ttV/ 2 
243 l 



+ 
+ 
+ 



3672 - 288vr 2 - 1296C(3) + (1944 - 324vr 2 )L + 108L 2 + 36L ; 



324 - 576^ 



'1728 - 216tt 2 )L + 324L 2 + 36L C 



4ni t 
"~8T i 



1296 - 12tt 2 + 1776L - 2052L 2 
144 - 6vr 2 + 18L + 18L 2 



+ 



-54 - 6vr 2 + 108L + 18L 2 



+ [116-96L] z 3 } + C{z^L 



(16) 



where L = logz. This function is renormalization scheme independent. All scheme depen- 
dent terms in the two-loop matrix elements of 2 and 2 undergo GIM cancellation and 
only affect the 0(y* b V ub ) part of eqn. (0). 

In the interesting range 0.005 < z < 0.1, the function f(z) is approximated (to 15% 
accuracy) by 



/(*) 



■19,2(1 - 2z) - 3^1og 2 z. 



;i7) 



2 Two-loop diagrams with 0\ insertions vanish due to their color structure. One-loop c — > 1*7 diagrams 
with 0\ or 0\ insertions vanish for on-shcll photons. 
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Using the s-quark mass m s {fi = m c ) = (0.17 ± 0.03) GeV || (but keeping m c = 1.5 GeV 
fixed), we find 

\A\ = \V; S V US \^^- (0.86 ± 0.19) = (4.7 ± 1.0) x 10" 3 . (18) 

All 

Thus, the two-loop amplitude is more than two orders of magnitude larger than the lead- 
ing logarithmic result in eqn. (|i~3|) and four orders of magnitude larger than the one-loop 
contribution in eqn. (|2|). The three contributions to the coefficient A are summarized in 
table 1. 

\AA lloop \ ~2 x 10- Y 



One-loop electroweak diagrams 



Leading logarithmic approximation 



Dominant two-loop diagrams 



\AA LLA \ = (8 ±3) x 10- 



\A\ = (4.7 ± 1.0) x 10" 



Table 1. Summary of the contributions to the coefficient A. 

The form of eqn. ( |I4D ensures us that no furher enhancement of the perturbative am- 
plitude is expected at even higher orders. Any AS = charm decay must be Cabibbo- 
suppressed. The suppression by m s /m c must also remain. The latter suppression is actually 
rather mild, as one can see from eqns. (|H)D and (|2"0"D below. We have thus exhausted the 



possibility of finding large contributions by considering higher orders in perturbation theory. 

4. The dominant two-loop contribution to the c — > wy amplitude is suppressed by (m s /m c ) 2 . 
Since m s is smaller (though not by much) than the A parameter of HQET, the nonperturba- 
tive contribution to a D-hadron decay may be equally (or more) important than the c — > wy 
contribution. 

Some part of the nonperturbative contribution may be taken into accont by replacing 
current quark masses in the r.h.s. of eqn. ([T4]) by constituent ones. Since the ratio of 
constituent masses is closer to unity than the ratio of current masses, one might worry that 
nonperturbative contributions may tend to bring back the GIM cancellations between s and 
d quarks. It is easy to explicitly check that this is not the case. For the constituent masses 
m c p = 0.3 GeV and m c s on = 0.45 GeV one obtains 

f[{mf n /m c ) 2 } - f[(m c d on /m c ) 2 ] ~ -0.68 - 0.48 i, (19) 

while for the current masses one has 

f[(m s /m c ) 2 ] - f[(m d /m c ) 2 ] ~ -0.24 - 0.68 i. (20) 
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We see that the resulting amplitude does not decrease (it is even larger) and conclude that 
the enhancement effect is rather robust. 

Ignoring the nonperturbative terms and normalizing the c — > uy rate to the semileptonic 
decay rate (analogously to what one usually does in the b — > sy case), we would obtain the 
following c — > uy contribution to the branching ratio of AS = weak radiative D-meson 
decay: 

ABR[ D X 7 l = ^ER^l BR{ D ^ Xev\. (21) 

7T \V cd \ 2 

This amounts to roughly 5 x 10~ 8 for D + and to 2.5 x 10~ 8 for D°. Although it is dramatically 
increased over the previous results, it remains very small. Compared to c —>■ uddy and c — ► 
ussy, the c — ► uy channel is suppressed (in branching ratio) by a factor of order \A/V c d\ 2 ~ 
5 x 10~ 4 . We conclude that there is little chance to extract the c — ► wy contribution even 
with use of the photon energy spectra. 

The obtained branching ratio is roughly confirmed when one looks at the exclusive decays 
D — > pry. Assuming its BR to be down by about a factor of ten from the inclusive one, the 
contribution of c — > wy amounts to an exclusive BR of about 5 x while the four- 
Fermi processes would give 10 -5 . This is in satisfactory agreement with a recent sum rule 
calculation of D — > pry. 

Since the c — > wy contribution in eqn. ( P*[j ) turns out to be suppressed by both (a s (m c ) /n) 2 
and powers of m s /m c , it is by no means the largest even among subdominant contributions. 
Decay modes with more gluons may have larger BR. For instance, the c — > wy gluon 
contribution to BR[ D A7] is proportional only to linear a s (m c )/TT. It can be written 
in the same form as eqn. ( pT|) but with \A\ replaced by 1.1 x 10 -2 . Consequently, it is a 
factor of 5 larger in the rate than the c — > wy contribution. The quoted numerical value 
of the c — > wy gluon contribution has been obtained with use of analogous b — > sy gluon 
results from ref. |[ITJ where appropriate replacements of electric charges had to be made. 



5. To conclude, we have pointed out an interesting enhancement of the (formally) higher 
order contributions in the c — ► wy decay. Nevertheless, the corresponding contribution to 
the inclusive weak radiative charmed hadron decay most probably remains screened by other 
channels. 

6. Note added: The present version of this paper differs from the original one by correcting 
a mistake in the negligible one-loop contribution in eqn. (fj). It was noticed by Q. Ho-Kim 



7 



and X.Y. Pham in 
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